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ABSTRACT 

We investigate the AdS / CFT correspondence for higher-derivative gravity sys- 
tems and develop a formahsm in which the generating functional of the bound- 
ary field theory is given as a functional that depends only on the boundary 
values of bulk fields. We also derive a Hamilton- Jacob i-like equation that 
uniquely determines the generating functional, and give an algorithm calculat- 
ing the Weyl anomaly. Using the expected duality between a higher-derivative 
gravity system and N = 2 superconformal field theory in four dimensions, we 
demonstrate that the resulting Weyl anomaly is consistent with the field the- 
oretic anomaly. 
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1 Introduction 



Over the past few years, many attempts have been made to check the AdS/CFT corre- 
spondence (For a review, see Ref. As an example, it is shown in Ref. 0] 
that the spectrum of chiral operators of A/" = 4 super Yang-Mills in four dimensions coin- 
cides with that of the Kaluza-Klein modes of type IIB supergravity on AdS^ x 5*^. Also, 
the computation of anomalies via bulk gravity has been shown to exactly reproduce the 
results of the super Yang-Mills theory However, this matching of the anomalies 
is valid only in the regime where N-^oo, X = gYM^^^j since the analysis is based on a 
classical supergravity computation. At present, it remains an important issue to test the 
duality beyond this regime. 

There have been several attempts to confirm the validity of the duality beyond the 
classical gravity approximation 0] HOI [|I31 • Among these, Ref. [§ treats N =2 G = 
USp{N) superconformal field theory (SCFT) in four dimensions. This SOFT can be 
realized on the world volume of D3-branes situated inside eight D7-branes coincident 
with an 07~ brane, and is known [|T^ to be dual to type IIB string on AdS^ x S^/Z2- 
The authors of Ref. showed that this duality reproduces the correction to the 
U{l)n chiral anomaly correctly. In Refs. and |jlO[, the correction to the Weyl 
anomaly of the SCFT is computed using a higher-derivative gravity theory in which a 
curvature square term is added. 

However, higher-derivative gravity theories^ exhibit some features in the AdS/CFT 
correspondence that differ from those in Einstein gravity. To see this, we first recall 
that the equation of motion for Einstein gravity is a second-order differential equation in 
time r. Thus, a classical solution can be totally specified by prescribing the value at the 
boundary if we further impose the regular behavior of the solution inside the bulk P], 
and the boundary value can be identified with an external field coupled to an operator 
in the dual CFT The situation changes drastically if we consider higher-derivative 

theories. In fact, a higher- derivative system with Lagrangian density C{g,g, ■ ■ ■ ,g^^~^^^), 
where gij is the metric and ■ = d/dr, generically gives an equation of motion that is a 
differential equation of 2(A+l)-order in r. We then would need (A+1) boundary conditions 
for each field to specify a classical solution, even if we require its regular behavior inside 
^For a review of higher derivative gravity, see, e.g., Ref. fl^. 
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the bulk. 

The main aim of the present paper is to formulate higher-derivative gravity systems 
in accordance with the holographic principle. In this paper, we say that the holographic 
principle holds when the following two conditions are satisfied: (1) the classical solution of 
a higher-derivative system is specified uniquely by the boundary value of each bulk field, 
and (2) the bulk geometry becomes AdS-like near the boundary. In order to satisfy the 
first condition, we first note that the system C{g^ 9, - ■ ■ , g^^~^^^) can be transformed into a 
Hamilton system with (A^+1) pairs of canonical variables [g, Q'^), {p, Pa) = 1, ■ ■ ■ ? ^) 
by defining Q"^ = d°'gij/dr"'. (See the next section for details.) Thus, by setting boundary 
conditions that are of the Dirichlet type for g and the Neumann type for Q", the classical 
solution of this system can be specified only by the boundary value of g. Note also that 
the classical action of this system, which is obtained by plugging this solution into the 
action, becomes a functional of these boundary values of bulk fields. The second condition 
ensures the existence of a UV fixed point of the dual theory at the boundary, and such 
a fixed point enables us to take the continuum limit We see below that appropriate 
boundary terms need to be added to the bulk action in order for the bulk metric to exhibit 
such asymptotic behavior when higher-derivative terms exist. 

For a systematic treatment of these issues, we employ the Hamilton- Jacobi formula- 
tion, as introduced by de Boer, Verlinde and Verlinde |TH] to investigate the holographic 



RG structure of Einstein gravity. (See Refs. [|T6|-p4| for more details of the holographic 
RG.) This formulation is further elaborated in Refs. [^-[0- In particular, a systematic 
prescription for calculating the Weyl anomaly in arbitrary dimensions is developed in Ref. 
p6| . In this paper, we show that the Hamilton- Jacobi equation is quite a useful tool also 
to study the holographic RG structure in higher-derivative systems. Actually, we can 
derive a Hamilton- Jacobi-like equation that determines the classical action in accordance 
with the holographic principle. That is, the classical action can be solved as a functional 
of a boundary value for each bulk field. As a check of our formulation, we compute 
corrections to the Weyl anomaly of the Af = 2 SCFT by solving the Hamilton- Jacobi-like 
equation. In the course of this analysis, we find that the prescription developed in Ref. 
p6| is again helpful. We show that our result can reproduce that of Refs. and |[T0[| . 

The organization of this paper is as follows. In §2, we formulate the Hamilton- Jacobi 
equation for a higher- derivative system with emphasis on applications to the AdS/CFT 
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correspondence. In §3, we apply the formulation to higher- derivative gravity and derive 
an equation that determines the classical action. In §4, we solve the equation following 
the prescription given in Ref. |2^, and demonstrate how to calculate the Weyl anomaly. 



We show that the resulting Weyl anomaly correctly reproduces that given in Refs. [Q] and 
10| . Section 6 is devoted to a conclusion. There, a comment is given on the holographic 



RG structure in higher- derivative gravity systems. Some useful results are summarized in 
the appendices. 

2 Hamilton-Jacobi equation for a higher- derivative 
Lagrangian 

In this section, we give a prescription for determining the classical action when higher- 
derivative terms are added. We start our discussion for a system of point particles with 
the action 

5[g(r)] = / drL(g,g,--- (g(")(r) = rf"g(r)/rfr") . (2.1) 

Jt' 

The extension of our argument to gravitational systems is straightforward and will be 
carried out in the next section.0 

The action (|2.1| ) can be rewritten into the first-order form in the following way. We 
first introduce the Lagrange multipliers p, Pi, - ■ ■ , Pn-i, so that q,Q^ = q, ■ ■ ■ , Q'^ = q^^^ 
can be regarded as independent canonical variables: 

= P{q - Q') + Pi{Q'-Q') + --- + PN-iiQ""-' - Q"") 

+ L{q,Q\--- ,Q'',Q''). (2.2) 

We then carry out a Legendre transformation from {Q^ , Q'^) to (Q^, Pn) through 

dL 



N 



(g,g\---,Q^g^) . (2.3) 



We here assume that this equation can be solved with respect to (= / (g, Q^, ■ ■ ■ , Q^; Pjy) ) , 
^See also Ref. ||3^ , where higher-derivative systems are discussed from the viewpoint of string theories. 
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and thus obtain the following action that is equivalent to (|2.1|) classically: 



Jt' 



dr 



N 



pq + Y,P-Q"-H{q,Q''\ p. Pa 



a=l 



(2.4) 



where Q is now the time-derivative of the independent variable Q , and the Hamiltonian 



is given by 



H{q,Q'^;p,Pa) = pQ^ + P,Q^ + ... + P^_,Q^ + P^f{q,Q-;P^] 
-L{q,Q\--- ,Q'^J{q,Q^;P^)). 



(2.5) 



The variation of the action ( |2.4| ) is given by 



6S 



and thus the equation of motion consists of the usual Hamilton equations, 



• dH ■ OH . 

dp dPa 



dH 



, Pa 



dH 



dq ' " dQ''' 

and the following constraint, which must hold at the boundary, r = t and r = t': 



pSq + Y.PaSQ" = {r = t,t') 



(2.6) 



(2.7) 



(2.8) 



The latter requirement, ( |2.8| ), can be satisfied when we use either Dirichlet boundary 
conditions. 



Dirichlet : 5q = , SQ" = {r = t, t') 
or Neumann boundary conditions, 

Neumann : p = , = {r = t,t') 



(2.9) 



(2.10) 



for each variable q and (a = 1, ■ ■ ■ , A^). If, for example, we take the classical solu- 
tion {q,Q°-,p, Pa) that satisfies the Dirichlet boundary conditions for all {q-iQ"") with the 
specified boundary values as 



q{r = t) = q, Q''{r = t) = Q\ and q{r = t') = g', Q"(r = t') = Q'" , (2.11) 
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then after plugging the solution into the action, we obtain the classical action that is a 
function of these boundary values, 



S{t, q, t', q', Q'^) = S [q{r), Q^ir); p{r), P,(r)] . (2.12) 

However, as we discussed in the Introduction, this classical action is not of great interest 
to us in the context of the AdS/CFT correspondence, since the holographic principle 
requires that the bulk be specified by only the values q and q' at the boundary. This leads 
us to use mixed boundary conditions: 

Sq = Pa = {r = t,t'). (2.13) 

That is, we impose Dirichlet boundary conditions for q and Neumann boundary conditions 
for In this case, the classical action (to be called the reduced classical action) becomes 
a function only of the boundary values q and q': 

S = Sit,q;t\q'). (2.14) 

A renormalization group interpretation of this condition is discussed briefly in the con- 
cluding section, and will be discussed in detail in a forthcoming paper |^ . 



Now we derive a Hamilton- Jacobi-like equation that determines the reduced classical 
action ( p.l4|) . This can be derived in two ways, and we start with the more complicated 



way, since this gives us a deeper understanding of the mathematical structure. To this end, 
we first change the polarization of the system by performing the canonical transformation^ 

S = S-[ dF, (2.15) 
Jt' 

with the generating function 

F = 5^P,Q\ (2.16) 



■^The following procedure corresponds to a change of representation from the Q-basis to the P-basis 
in the WKB approximation: 
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Although the Hamilton equation does not change under this transformation, the boundary 
conditions at r = t and r = t' become 

p5q-Y,Q"^Pa = ^ {r = t,t'). (2.17) 

a 

These boundary conditions can be satisfied by imposing the Dirichlet boundary conditions 
for both q and Pa'- 

q{r = t) = q, Pair = t)=Pa, and q{r = t') = q\ Pair = t') = P;^. (2.18) 

Substituting this solution into S, we obtain a new classical action that is a function of 
these boundary values, 

S{t,q,Pa, t',q',P:^) = S [q{r),Q-{ry, p(r),P,(r)] . (2.19) 

By taking the variation of S and using the equation of motion, we can easily show that 
the new classical action S obeys the Hamilton- Jacobi equation: 

dS / dS dS ^ 

m=-'' ['^-dPa^ 

9S ,J , dS OS \ 

The reduced classical action S{t, q; t', q') is then obtained by setting Pa = in S: 

S (t, q- 1', q') = S (t, g, = 0; t', q', P^ = 0) . (2.21) 

Note that the generating function F vanishes at the boundary when we set Pa = 0. 

Here we briefiy describe how the Hamilton- Jacobi equation ( p. 201) is solved. For sim- 
plicity, we consider the case = 1 and focus only on the upper boundary at r = t. 
Motivated by the gravitational system considered in the next section, we assume that the 
Lagrangian takes the form 

Hq, q) = Loiq, q) + cLi(g, q, q), (2.22) 

where 



Lo{q,q) = -mij{q)q'q^ - V(q), 



Li{q, q, q) = ^n,j{q)q'q^ - Ai{q, q)q' - (f){q, q), (2.23) 



with 



Hq, q) = &mq'q'q' + 0!?(g)^Y + 0^°^?)- 



(2.24) 



We further assume that the determinants of the matrices rnij{q) and nij{q) have the same 
signature.0 Following the procedure discussed above, this Lagrangian can be rewritten 
into the first-order form 



L = pq + PQ~Hiq,Q;p,P) 



(2.25) 



with the Hamiltonian 



H{q,Q-p,P) = p,Q^ - -m,j{q)Q^Q^ + V{q) 



+ Y^n^'ii) ( + ''Mq: Q) )[Pj + cMq, Q) ) + c</)(g, Q), (2.26) 



where n*-' = {n.ij) ^. The Hamilton- Jacobi equation ( p.20|) is solved as a double expansion 
with respect to c and P by assuming that the classical action takes the form 

1 



S{t,q,P) 



S^,/2{t, q, P) + 5o(t, g, P) + v^Si/2(t, g, P) + c Si{t, g, P) 



+0{c 



3/2 ^ 



^2.27) 



After some simple algebra, the coefficients are found to be 
1 



= -u'^iq)P.P, + OiP'), 
So = So{t,q)~P,d'So + 0{P' 



5*1/2 = PiU'^{q)njk{q) 
+C(P2). 



rL d'So d^So + d'V{q) + n'\q)Ai ( g. 



dq 



(2.28) 



Here, 



d_ 

(9g* 



(9* = m}^dj 



(2.29) 



and r*^ is the affine connection defined by m^. Also u^^ is defined by the relation 



u"'{q)u^\q)mki{q) = n'^{q). 



^2.30) 



*In fact, it is easy to see that this is the case in the higher-derivative gravity system considered below. 
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Furthermore, So{t,q) = So{t,q, P = 0) and Si{t,q) = Si(t,q, P = 0) satisfy the equations 
dSo 



dt 



1 , .dSo dSo , , 



mij{q) 



dSidSo 
dq^ dqi 



-n,,{q) (rL d'S, d'So + dW{q)) [V^^ d^S, d^S, + &V{q)) 



-AAq 



{Tlid'Sod^So + d'V{q))+Mq 



dSn 



(2.31) 



which can be expressed as a Hamilton- Jacobi-hke equation for the reduced classical action 
S{t,q) = So{t,q) + cSi{t,q) + 0{c^y. 

dS 



where 



H{q,p) = -m'^{q)piPj + V{q) 



Pi 



[2.32] 



+ c 



- AM.P) (n^pV + d'V{q)) + <P{q,p) 



(2.33) 



It is important to note that H is not the Hamiltonian. In fact, the Hamilton equation for 
H does not coincide with that obtained from ( p.26|) . 



In solving the full Hamilton- Jacobi equation ( p.20| ) for S{t,q,P), we imposed the 
condition that everything becomes regular around c = when we set P = 0. This is 
because in most interesting cases (like those of the gravity systems we discuss in the 
following sections) the higher-derivative term is regarded as a perturbation, so that the 
reduced classical action must have a finite limit for c^O. Once such a regularity condition 
is imposed, we have an alternative way to derive this pseudo-Hamiltonian H with greater 
ease. In fact, for any Lagrangian of the form 



L{<t,cf,ct) = Lo{<t,ct) + cL,{(t,cf,if) 



(2.34) 



one can prove the following theorem, assuming that the classical solution can be expanded 
around c = 0:| 

^As long as we think of Li as a perturbation, any classical solution can be expanded as 

q{r) = qo{r) + cqi{r) ^ 0{c^) . 
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Theorem 

Let Ho{q,p) be the Hamiltonian corresponding to LQ{q,q). Then the reduced classical 
action S{t,q; t',q') = So{t,q; t' ,q') + c Si{t,q] t',g') + 0(c^) satisfies the following equation 
up to 0{c^): 

dS Tr/ N dS ^ dS ~, , n , dS 

-^ = H{q,p), P. = ^, and + — = H{q\p% K = (2-35) 

where 



H{q,p) = Ho{q,p) - cLi{q, fi{q,p), f2{q,p)) 
fliq,p) ^ {Ho,q'} 



r,{q,p) ^ {Ho,{Ho,q'}} 



dpi ' 



dpidq^ dpj dpidpj dq- 



jpr \r( ^x-dPdG dGdF\ 
{F(g,p),G(g,p)} = -----J (2.36) 

A proof of this theorem is given in Appendix A. It can easily be confirmed that this 
correctly reproduces (p.32|) and (|2.33|) for the Lagrangian given in (p.22|)-(^.24|). 



3 Application to higher- derivative gravity 

In this section, following the prescription developed in the previous section, we derive 
an equation that determines the reduced classical action for a higher- derivative gravity 
system. 

We first recall the holographic description of RG flows in the dual boundary field 

theory. We parametrize the bulk metric with the Euclidean ADM decomposition. (For 

Here go is the classical solution for Lq, and qi is obtained by solving a second-order differential equation. 
Note that we can, in particular, enforce the boundary conditions 

qo{r = t)=q, qi{r = t) = and qo{r = t')=q', qi{r = t') = 0. 

In this case, due to the equation of motion for qoir) , the classical action is simply given by 

S{q,t;q',t') = dr[Lo{qo, %) + c Li{qQ, +0{c^). 

This corresponds to the classical action considered in Ref. ||l0|] . 
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more details of the ADM decomposition, see Appendix B.) We then have 

= N{x,rYdr'^ + gij{x,r){dx' + \\x,r)dr^{dx^ + \^{x,r)dr^. (3.1) 

Here X'^ = (x*,r), with i,j = 1,2, ■ ■ ■ ,d, and N and A* are the lapse and the shift 
function, respectively. The signature of the metric Qij is taken to be (+■■■+). By 
assuming that the geometry becomes AdS-like in the limit r-^— oo, the Euclidean time r 
is identified with the RG parameter of the d-dimensional boundary theory, and the time 
evolution of other bulk fields (such as scalars) is interpreted as an RG flow of the coupling 
constants with a UV fixed point at the boundary. To avoid a singularity of the metric g^j 
at r = —oo, we restrict the region of r such that rg < r < oo [00||3^. This corresponds 
to the introduction of a UV cut-off to the boundary field theory. In the following, we 
consider a (c? + l)-dimensional manifold M^+i = {{x^,r)} that has a topology given by 
Md+ir^{R^ U oo) X R+, with ro<r<oo. 

We consider classical gravity on M^^i with the action 

S = SB + St. (3.2) 
Here Sb is the bulk action given by 

Sb = I d'^'-'X^CB, (3.3) 
Lb = 2A- R-aR'-bR^-cRl^^, (3.4) 

where a, b and c are some given constants. Sb contains boundary terms defined on the 
boundary = dMd+i at r = tq. The form of Sh can be determined by requiring that it 
is invariant under the diffeomorphism 

Xi" X'^" = /^(X), (3.5) 

with the condition 

r {r = ro, x) = ro. (3.6) 



Equation (|3.6|) implies that the diffeomorphism does not change the location of the bound- 
ary. It is then easy to verify that Sb takes the form (for details see Appendix C) 

d'x^B, (3.7) 



d 
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with 



(3. 



where Kij is the extrinsic curvature of Y^d given by 



2N ^^'^ 



(3.9) 



and K = g^^Kij. Vj and Rijki are, respectively, the covariant derivative and the Riemann 
tensor defined by gij. The first term in B ensures that the Dirichlet boundary conditions 
can be imposed in the Einstein theory ||3^ and also plays an important role in the context 
of the AdS/CFT correspondence |3^. We argue below that the coefficients xi, ■ ■ ■ ,0:5 must 
obey some relations so that the holography holds even for higher-derivative gravity.^ 

The action (|3.2| ) is expressed in terms of the ADM parametrization as 

d 



whereQ 



Md+i 



ro 



dr / d'^x^ Co{g, K- N, A) + Ci{g, K, K; N, A) 



(3.10) 



-Co = 2A-R + ^-K' 



(3.11) 



N 



-aR^ - bR 



6a + 2xi)Kf. + (2a - xi)K^ 



R 



+ 



2(26 + 4c - X2)iK^)ij + (26 + 2xi - X2)KK, 



R'- 



+ 2{6c + X2)K,kK,iR'"'' 



- 2(26 + c - 3x5)Kfj + (46 + 4x4 - X5)KKfj 

6x3 



- (9a + 6 + 2c - 2x4) {Kfj) + (6a - 6 

- (a + X3)K^ 



Xi)K^Kl 



- (46 + 2xi - X2)K,jV'V^K + 2(6 - 4c + X2)K,jV^VkK 



ki 



^See, e.g., Refs. and |Q for another discussion of boundary terms in higher-derivative gravity. 
^We here use the fohowing abbreviated notation: AT^ = Ki^K}!^ ■ ■ ■ K^, {K%j = K^kK^ . 
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+ (8c + X2)K,jV^K'^ + 2(6 + xi)KV^K 
- [(4a + b)g'^g''' + {b + Ac)g''g''] U^L^i 

|(4a - Xi)R + (12a + 26 - Xi)Kli - (4a + ?>x^)K'^^g'^ 
+ (26 - X2)R'^ + (46 + 8c - 'ix^){K^ f - 2(6 + Xi)KK'^ 



+ 



(3.12) 



with 



(3.13) 
(3.14) 



By regarding gij and Kij as independent canonical variables,^ the action ( |3.10| ) can be 
further rewritten into the first-order form 



/•oo I* 

S= dr d'^x^ [n'^ {g,, - - V,A, - V, A,) + £o + ^i] 

J rn J 



ro 
oo 



dr J d^Xy/g yji'^gij + P'^ K,j - n{g, K; vr, P; N, X) 

Here the Hamiltonian density 7t can be evaluated as 

n = 7i'^ {2NKij + ViXj + VjXi) + P'^kij -Co-Ci 
= Nn{g,K-7r,P) + X,V\g,K-7i,P), 



(3.15) 



(3.16) 



with 



n{g,K;7r,P) = 2n'^ K, 



-P.2. + 

.7.1 I 



4a + 6 



4(6 + 4c)^*^ ' 4(6 + 4c)(4rfa + (d+l)6 + 4c)' 
- WiVjP'' + [Ai R'^ + A2 {K^f + A3 KK'^] Pij 

+ [AiR + A5Kf^+AeK^]P 

-2A + R- Kl + K"^ 



+ 5i R^ + B2Rl + B^R^ 



ijkl 



®The correspondences between the variables in §2 are as follows: q ^ gij, p ^ *5 ^ ^ 
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+ (Ci Kl + C2 K') R + [C3 iK%^ + C,KK,,] 



-kl 



+ KijV^K'^ + irV^ir, (3. 17) 



TT, P) = - 2Vy' + PkN'K'' - 2V'{K''Pjk). (3.18) 

The coefficients Ai, ■ ■ ■ , i?4 are not important in the following discussion, and are listed 
in Appendix D. The classical equivalence between the two actions (|3.10|) and ( 3.15 ) can 



be easily established by noting that the latter gives the following equation of motion for 

P'^ = - 2 ((4a + b) g'^g'' + (6 + 4c) (7 V) L^i 

+ [(4a - xi)R + (12a + 2b - Xi)Kli - (4a + 3x3)i^^] g'^ 

+ {2b - X2)R'^ + (46 + 8c - ?>x^){K^y^ - 2{b + X4) KK'^ . (3.19) 

This correctly reproduces the original action ( |3.10| ) when substituted into ( p.l5| ). 

Following the prescription given in §2, we now make a canonical transformation that 
changes the polarization of S from {gij, Kij) to {gij, P*-'): 

S^S- I d'^'X^i^K^.P^^) 

= ^ jd'^x^ (n'^g,, - K,.jP'^ -NH- X.V') , (3.20) 



with 



n{g, K- TT, P) = n{g, K- vr, P) + KK,,P'^, 
V\g, K; vr, P) = V\g, K; vr, P) - V^K.kP'") 

= -2Vjn'^ - V'P^'' Kjk - 2V''{K'^Pjk), (3.21) 

where we have used the relation 

dry^=^{NK + V'X,). (3.22) 

Since N and Aj are the Lagrange multipliers, we obtain the Hamiltonian and momentum 
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constraints 



J_5S_ 
\f9 ^\ 



n{g,K- 71, P) = 0, 
V\g,K; 71, P) = 0. 



(3.23) 
(3.24) 



We now let Qij and P*-^ represent the solution to the equation of motion for S that 
obeys the boundary conditions 



gij (x, r = ro) = gij{x) , P'^{x,r = ro) = P'^ (x) . 



(3.25) 



We also require that the solution be regular or be set to some specific value inside the 
bulk (r— s>oo), and assume that the above boundary condition is sufficient to specify the 
classical solution completely 0. Plugging the solution into S, we obtain the classical 
action S[g{x), P{x)], which satisfies the following Hamilton- Jacobi equation:]^ 

1 6S ,, 1 6S 



-71 



+K, 



ij , 



H{g,K-7,,P) 
V\g,K-7i,P) 



0: 
0, 



(3.26) 

(3.27) 
(3.28) 



Since the Hamiltonian density is a linear combination of the constraints, the classical 
action S does not depend on the coordinate of the lower boundary: 



This implies that the reduced classical action 



S[g{x)]^S[g{x),P{x)=Q] 







(3.29) 



(3.30) 



^The last equation demonstrates the invariance of S under a c?-dimensional difTeomorphism, 



= - / d'^x^e.V 



(V.e, + V,eO^^ + {-due' P''^ - dke^ P^'^ + e'^ dkP'') 



with e'(x) an arbitrary function. This also demonstrates the invariance of the reduced classical action, 

f S S 

= / d'^x(V,ej+V,e,) 

for arbitrary e*(x). 
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is also independent of tq: 



(3.31) 



The Hamiltonian and the momentum constraints ( |3.27| ) and ( |3.28| ) can be translated 
into equations for the reduced classical action, as we sketched for point-particle systems in 
Eqs. ( p.27| )-( p.33| ). However, the resulting equation can be derived most easily by using 
the Theorem, ( p.35|) and (|2.36|) , as follows: We first rewrite the Lagrangian density of 
zero-th order, Co, into the first-order form 



where the zero-th order Hamiltonian density TYq is given by 

1 



(3.32) 



Uoig, it; N,X) = N \^7rf^ - j-^tt' - 2A + i?J - 2A, V.vr*^ . 
Then by using the Theorem, the pseudo-Hamiltonian density is given by 
nig, tt; N, A) = Hoig, vr; N, A) - ^g, K\g, vr), K\g, vr); N, A) . 



(3.33) 



(3.34) 



Here K^-{g, vr) is obtained by replacing gij{x) in ( |3.13[ ) with | J d'^y^l-Co{y), (^^(x)}, and 
it is calculated to be 



1 



(3.35) 



On the other hand, K^j = {/ d'^yy/gT-Co{y), Kfj^ is found to be equivalent to replacing 
Lij in Ci by 

1 



7-0 ^ 

2{d-iy 



2{d - 1)A + {d- 1)R +{d- - 37r' 



9ij 



+ Rij + 2{7l'^)ij - ^ _ ^ TTTlij . 



(3.36) 



Using Eqs. ( p.31| )-( p.36| ), we obtain the following Hamilton- Jacobi-like equation for the 
reduced classical action: 



= j d''x^n{g{x),7c{x);N,X') 

d'^x^ \Nn{g{x),n{x)) + X' V^{g{x),^l{x)) 



(3.37) 
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where0 



with 



n{g,7r) = 



d- 1 



-1 5S 



TT^ -2A + R 



TT 



2 



+ 7i + 72 Ai? + 73 + 74 Rl + 75 i?^^ 



«! = 2c, ^2 = 
1 

~4((i- 1)2 



2x, 



(rf-1)' 

4a + (d^ -3d + 4)6 + 4(rf - 2)(2d - 3)c 



04 



"5 



-2{d-l){dx4 + 3x5) 

, } r-4a - id^ - 3c/ + 4)6 - 4(2c/2 - 5rf + 4)c 
2(a — 1)'^ 

- 3rfx3 + (2c/2 - 7rf + 2)x4 - 3(2rf - 1)^5] , 
[4a + {d^ - 3c/ + 4)6 + \{2d^ - 5d + 4)c 

+ 2(3c/ - 4)x3 - 2(c/2 - 6c/ + 6)^4 + 2(5c/ - 6)x5] 



(3.38) 



(3.39) 
(3.40) 



(3.41) 



""^^We have ignored those terms in Ti, that contain the covariant derivative V. This is justified when 
we consider the holographic Weyl anomaly in four dimensions. Actually, it turns out that they give only 
total derivative terms in the Weyl anomaly. 
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/34 



1 



Ada - d{d - 3)6 - A{d - 2)c - (rf - l)(rfx4 + 8x5 

- 4da + d{d - 3)6 + A{d - 2)c 

- 3da;3 + {d^ -2d- 2)x4 + 3{d - 2)x5 
4a + {d^ -3d + 4)6 - 4(3rf - 4)c 

- (d — l)((ixi + X2 — (c? — 2)x4 + 3x5) 

- 4a - -3d + 4)6 + 4(d - 2)c 

- (d - l)(rf - 4)a;i - 3(rf - l)x2 + 3{d - 2)x3 

- {d^ -8d+ 10)x4 + 3(3rf - 4)0:5 

/?5 = 16c + 3X5, A 



{d-iy 

1 



1 

2(rf- 1) 



1 

2{d-l] 



2{xi + 2x2 — Xi — 3x5) 



d-l 



Pi 



12c -2x2, 



(3.42) 



d 



71 

72 
73 

74 



(rf-l)2 
1 

1 

4(d- 1)2 

4c + X2, 



4da + (rf + 1)6 + 4c , 

Ada - d{d - 3)6 - A{d - 2)c - {d - l){dxi + X2) 

4a + {d^ -3d + 4)6 - A{3d - A)c + 2{d - l){{d - 2)xi - X2) 
75 = c. 



(3.43) 



Since the classical action S[g{x), P{x)] is independent of the choice of and A* (and, 
thus, so is S[g{x)]), from Eqs. ( p.37|) -( ^.40|) we finally obtain the following equation that 
determines the reduced classical action: 

-1 5S 



n{g,j{x),7r'^{x)) = 0, V,{g,j{x),7r'^{x)) = 0, 7r*^(x) 



(3.44) 



We conclude this section by making a few comments on the possible form of the 
boundary action Sb and the cosmological constant A. As discussed above, in order that 
the boundary field theory has a continuum limit, the geometry must be asymptotically 
AdS: 



ds' 



dr^ + e '^'^^'''yij{x)dx^dx^ for r ^ —00. 



(3.45) 



This should be consistent with our boundary condition P*-' = 0. By investigating the 
equation of motion derived from the action ( p.l5|) explicitly, it can easily be shown that 



18 



this compatibility gives rise to the relation 

xi = 4a, 
X2 = 2b, 



X3 + dxi + X5 = —-(^d(d + l)a + db + 2cj . 



(3.46) 



It can also be shown that the asymptotic behavior ( p.45| ) determines the cosmological 
constant A as 



^_ d{d-l) ^d{d-3) 



2/2 



21^ 



d{d+ l)a + db + 2c 



(3.47) 



4 Solution to the flow equation and the Weyl anomaly 



In this section, we solve the equation ( |3.44| ), using the derivative expansion that was 
developed in Ref. p6|. We then apply the result to computing the holographic Weyl 
anomaly of A/" = 2 superconformal field theory in four dimensions, which is dual to IIB 
supergravity on AdS^ x S^/Z2- 

We first note that the basic equation, (|3.44| ), can be rewritten as a fiow equation of 
the form 



with 



i^f {S, S} 



{S, S} + {S, S, S, S} — Cd, 



(4.1) 



6S 



d-1 



6S 



5S 



9ij 



P2 A Qij 



6S 



6S 



+ /34 -R gij 



6S 



+ k RijQkl 



6S 6S 



(55* SS 



SS 6S 



kl 



Sgik Sgji 



(4.2) 



{^)'{S,S,S,S} ^ 



6S 



Sgij 
+ "4 ( gki 



+ "2 gki 



6S\ f 5S 



^gkij \Sg 



"3 



6S 
Sgij 



SgkiJ \Sg 



+ "5 gij 



6S 
' Sgij 



(4.3) 
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(4.4) 



Following Refs. |0 and ||2^, we then assume that the reduced classical action S[g{x)\ 
takes the form 

1 , 1 



S[g{x)] 



5'ioc[fi'(a;)] + T[g{x)\ 



(4.5) 



where 2k^_,_^ is the ((i + l)-dimensional Newton constant. The functional T[g] is identified 
with the generating functional of the boundary field theory in the background metric 
gij{x), with any local sources set to zero, and 5'ioc[fi'] is the local counterterm in S[g]: 



Sioc[gix)] 



(fx \/ g{x) Cxoc{x) 
fd'x^/^) [^ioc(x)]^. 



(4.6) 



«)=0,2,4,-' 



Here we have arranged the sum over local terms according to the weight w which is 
defined additively from the following rule: 





weight 


gij{x), T[g] 







1 




2 


ST/6gij{x) 


d 



We then substitute ( f4.5|) into the flow equation ( [4.1| ) and rearrange the resulting equation 
with respect to the weight. The parts of weight and 2 give 

2A - 7iA2 



-R - 72A/? 



{ 'S'loc ) 


•S'loc} 


+ 







{ "S'loc ; 


•Sloe} 


+ 

2 



{*SloC) *Sioc, "S'loc, "S'loc} 
{*SloC) *Sloc, >Sloc, >Sloc} 



J 2 



These two equations determine [i2ioc]o and [£ioc]2 as 



locJo 



W. 



10CJ2 



(4.7) 
(4.8) 

(4.9) 



11^ 



d-2 {d-l){d-2)l 

3{d'^X3 + dx4^ + X5) 



4:d{d + l)a — Mb — 8c + d{d'^X3 + dx^ + X5) 
d{d+ l)a + db + 2c 



1 

+ T 



dxi + X2 + 



2{d 



(4.10) 
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where (|3.47|) has been used. It is worthwhile to note that W and $ can be written in 
terms of only a, b and c upon substituting into ( p.46| ): 

2{d-l) 4((i + 3) 



W 



/ 2 



d{d + l)a + db + 2c 



d{d - 5)a - 2b - 2c 



(4.11) 



d-2 {d-2)l 

For d>4:, the flow equation of weight 4 simply determines [-Ciocj^ in the local counterterm, 



as in the case of Einstein gravity (cf. Ref. @), while for d = A this gives an equation that 
characterizes the generating functional r[(yf(x)]: 



{'S'loc, r} 

1 



J 4 



{•S'loc; "S'loc, *Sloc5 r} 



J 4 
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{•S'loc, "S'loc} + {"S'loc, "S'loc, "S'loc, "S'loc} 
J 4 L 

+ 73i?' + 74i??, + 75i?-,M)- 



(4.12) 

From this, we can evaluate the trace of the stress tensor for the boundary fleld theory: 

(4.13) 



(T.)g = —9ij- 



In fact, using the values in ( [4.10|) , we can show that the trace is given by 



2/3 

^^'^^ " 2^ 



— 1 5a b c \ 2 5a b 3c \ 2 



I ^ p2 

2P ^^^^ 



(4.14) 



This correctly reproduces the result[3 obtained in Refs. and |T0|, where the Weyl 
anomaly was calculated by perturbatively solving the equation of motion near the bound- 
ary and by looking at the logarithmically divergent term, as in Ref. H. 



^^The authors of Refs. p and [|oj parametrized the cosmological constant A as 

d{d-l) 



A = -- 



2L2 



so that their L is related to our Z, the radius of asymptotic AdS, as 

(d-3) 



1 - 



(d- 1)L- 



■{d{d+l)a + db + 2c) 
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For the case oi M =2 superconformal USp{N) gauge theory in four dimensions, we 

VoK^VZa) (radius of ^VZa)^ 



choose such that 



1 

24 



2^2 



(4.15) 



where 2^"^ = {2'Ky g1 is the ten-dimensional Newton constant |3^, and the radius of 
/ Z2 could be set to {SngsNy/'^ p|. In this relation, we note the replacement — 2N 
as compared to the AdS^ x 5*^ case. This is because here we must quantize the RR 
5- form flux over S^/ Z2 instead of over [0. For the AdS^ radius /, we may also 
set / = {SirgsNY'^. Setting the values a = 6 = and c/2P = 1/32N + 0{1/N^), as 
determined in Ref. 0], we find that the Weyl anomaly ([4.14|) takes the form 



N 



2 r 



27r2 



-1 1 
24 ^ 48iV 



1 



8 32N 



1 



p2 

32jv 



+ 0{N^). (4.16) 



This is different from the field theoretical result 



i 19 



27r2 



-1 

24 



32A^ 



W + 



1 1 

+ 



8 16A^ 



p2 I ^ 



+ C»(Ar°). (4.17) 



As was pointed out in Ref. [|T0], the discrepancy could be accounted for by possible 
corrections to the radius / as well as to the five-dimensional Newton constant. In fact, if 
these corrections are 

Vol(5VZ2) (87r(7,Ar)5/4 



2/t2 



(4.18) 



then the field theoretical result is correctly reproduced for 3^ + 77 = 5/4. 



5 Conclusion 

In this paper, we investigated higher-derivative gravity systems in the context of the 
AdS/CFT correspondence. Although higher-derivative gravity requires more boundary 
conditions than Einstein gravity, we pointed out that by choosing the Neumann boundary 
conditions for higher- derivative modes, the classical action can be made such that it 
depends only on the boundary values of bulk fields. We further derived a Hamilton- 
Jacobi-like equation that determines such a classical action. Using this equation, we 
computed the 1/A^ correction to the Weyl anomaly of M =2 G = USp{N) superconformal 
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field theory in four dimensions on the basis of the holographic description in terms of 
type IIB string theory on AdS^ x / Z2 [l^. We found that the resulting Weyl anomaly 
correctly reproduces the holographic Weyl anomaly given in Refs. § and JlOl, and is 
consistent with the field theoretical result if we take into account the possible corrections 



discussed in Ref. ||T0|| . 

Finally, we comment on how our Neumann boundary condition P = can be inter- 
preted in the context of the holographic RG. To this end, we consider a toy model with 
the Lagrangian 



(5.1) 



whose first-order form reads 



L = pq + PQ-H{q,Q-p,P), 



(5.2) 



with 



l ,2^2 1,02 



1 



H{q, Q; p, P) = -^fiY - ^Q' + QP + fP'- 



By performing an almost diagonal canonical transformation. 



with 



Q 
p 

\p J 



( 



1 



1 



fll ^2 ;^'^3 m2'^4 

C13 C14 (1\ C12 

cM^Os crn^a^ cM'^ai cm'^a2 



Q' 
p' 

\P' ) 



m = ^ = /i (1 + 0[c)j 



2c 



i±4^4(l + 0(c)). 



(5.3) 



(5.4) 



(5.5) 



^«3 + 



- 1 - 2c/i2 ' ^ M2 4 1 - 2c/i2 ' 
the Lagrangian can be rewritten into the following form with normalized kinetic term: 



(5.6) 



L=p'q' + P'Q'-H'{q\p'; Q\P') 



(5.7) 
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where 

H'{q\ Q'; p', P') = ^p" + \P'' - \m\'' - \m'Q'\ (5.8) 

Since a bulk mode with mass M is coupled to a scaling operator with scaling dimension 
A = I (d + Vrf + 4M2) 011], the relation shows that the mode Q' is coupled to a 
highly irrelevant operator with large scaling dimension when 1. The essential point 
of this conclusion does not change even if the variable q corresponds to a bulk field with 
spin. 

Turning to higher-derivative gravity systems, the above example shows that Kij (~ Q ~ 
Q') is highly irrelevant in the dual CFT and is approximated well by assuming that it takes 
a constant value along the renormalized trajectory, as long as we consider the vicinity of 
the conformal fixed point. This is equivalent to demanding that the corresponding beta 
function vanishes along the renormalized trajectory. Since P*-', the conjugate momentum 
of Kij, can be regarded as the RG beta function of Kij, this leads to our requirement, 
pij _ g fjj-^g holographic RG structure in higher-derivative systems will be explored in 



more detail in a subsequent paper 34 
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A Proof of Theorem 



In this appendix, we give a detailed proof of Theorem, ( p.35|) and (|2.36| ), for the action 



dr 



L,{q\ct) + cL,{q\q\it) 



(A.l) 



where i runs over some values. In the following discussion, we focus only on the upper 
boundary, for simplicity. 
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We first rewrite the zero-th order Lagrangian Lq into the first-order form by introduc- 
ing the conjugate momentum poi of as 



dr 



Poiq' - Ho{q,Po) + cLi{q,q,q) 



through the Legendre transformation from {q,q) to {q,Po) defined by 

POi 



dq 

From this, the equation of motion for poi and g* is given by 



q' 



POi 



dpi 



— + c 



dq' 



dLi d ( ^L^ \ d'^ f dLi 



(A.2) 

(A.3) 

(A.4) 
(A.5) 



(9g* dr \ dq' J dr'^ \ dq' 

Let q{r), po{r) be the solution to this equation of motion that satisfies the boundary 
condition 



q\r = t) = q' 



(A.6) 



Since this condition determines the classical trajectory uniquely [together with the lower 
boundary values q^{r = t') = q" that we have not written here explicitly], the boundary 
value of pq is completely specified by t and q: po{r = t) =Po(t, q). By plugging the classical 
solution into the action 5, the classical action is obtained as a function of the boundary 
value q' and t: 



S{t,q) = S[q{r),po{r)]. 



(A.7) 



In order to derive a differential equation that determines S{t, q), we then take the variation 
of Sit, q). Using ( |A.4| ) and (|A.5| ), this is easily evaluated to be 



5S 



poiq' - Ho{q,po) + cLi{q,q,q) 
dLi 



5t 

+ 5q\t) 



Poi + c 



dq^ 



iq,q,q) 



d ( dLi ^ 

[q, q) 



dr \ dq' 



r=t/ J 



where 



+ cSq\t) ^{q,q,q), 



dr'^ 



q' 
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(r = t) 



(A.8) 



(A.9) 



and Sq^{t) and Sq{t) are understood to be Sq^{r)\r=t and dSq\r)/dr\r=t, respectively. By 
expanding the classical solution g*(r) around r = t, we find that the variations Sq^{t) and 
6q{t) are given by 



Sq^ (t) = Sq' -q'St, (t) = Sq' - q' St. 



(A.IO) 



Here it is important to note that q can be written in terms of q and t, since the classical 
solution is determined uniquely by the boundary value q. Actually it can be shown that 



5q' 



dq^dpoi 



Sq^ + 
5q^ + 



O'^Hq [ dpoj_ 
dt 



6t + 



(A.ll) 



dqWpoi dpoiPoj 

where we have used ( |A.4j) as well as the fact that po = po(t, q). From these relations, the 
variation ( |A.8| ) is found to be 



5S=Pi Sq'-H{q,p) 5t, 



(A.12) 



with 



Pi 



Poi + c 
+ 



dLi 



riT \ rin'' 



+ 



dr \ dq 
d^Ho dp, 



r=t 



Ok 



dqi ydq'dpoj dpojdpok dq' 
H{q,p) = Ho{q,po) 



(A.13) 



+ c 



d ( dLi ^ 



Li{q, g, q) + q' -WT-{q, - 3- lFr(^' ^' 1) 



dq 

+ ^ (qi ^'^0 dpoj 



dr \ dq 



r=t 



(A. 14) 



dq' V dpoidpoj dt 

In order to compute H{q,p), we first note that the Hamilton equation appearing in ( [A .41 ) 
and ( |A.5|) gives the relation 



d^Ho dHo , d^Ho fdpojdHo , dp 



+ 



+ 



50fc 



dt 



dpoidq^ dpoj dpoidpoj \ dq'' dpok 
It is then easy to verify that H{q,p) takes the form 

H{q,p) = Ho{q,p) ~ cLi{q,q,q) + 0{c^). 

Here and g* in Li can be replaced by 



niq,p)^{Hoiq,p),q'} 



dHo 

dpi 



{q,p) 



(A.15) 



(A. 16) 



(A.17) 
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and 



n{q,p) ^ {Hoiq,p),{Hoiq,p),q'}} 



d^H dH d"^ H OH 



dpidq^ ' dpj ' dpidpj ' Sg^ 
respectively, up to C(c^). This completes the proof of ( p.35|) and ( |2.36| ). 

B ADM Decomposition 

In this appendix, we summarize the components of the Riemann tensor, Ricci tensor and 
scalar curvature written in terms of the ADM decomposition. 

In the ADM decomposition, the metric takes the form 

= N{x,rfdr'^ + gij{x,r)(^dx' + X'{x,r)dr^(^dx^ + X^{x,r)dr^. (B.l) 
Here we use the following basis instead of the coordinate basis (9^: 

en = ^{dr - X'di, ), = di. (B.2) 



(B.3) 



In this basis, the components of the metric are given by 

g(en,en) g(en,ej) \ _ I ^ 
giej,e^) g{ei,ej) J \^ gij 

For the purpose of computing the Riemann tensor in this basis, it is useful to start with 
the formula 

Ve,,Ve.] e^- V[e„e.]ep. (B.4) 

Each component can be calculated explicitly by using the equations 























= -^g'^'duNeu 















(B.5) 
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where K^j is the extrinsic curvature and F'j^. is the affine connection with respect to gij. 
We thus obtain 

Rijkl — Rijkl — KikKji + KiiKjk, 
Rnjkl — ^iKjk — VkKji, 

Rnjnl = - Ljl, (B.6) 

with 

Lij = 1 (^kij - X'^ VkK.j - Kkj - VjX' Kkj + V,V,-7V) . (B.7) 
The components of the Ricci tensor i?^,^ = R'^np,, — Rvfx are given by 

Rij ~ Rij ~ K^ij ~ Li^j, 

Rin = V'Kki - V,X, 

Rnn^^-g^'Lij, (B.8) 
and the scalar curvature is 

R^R + SKf.-K'' -2g'^Lij. (B.9) 

C Boundary Terms 

In this appendix, we supplement the discussion of the possible boundary terms given in 
§3. 

We first consider the infinitesimal transformation 

x' ^ x'' ^ x' + e'{x,r), r ^ r' ^ r + e{x,r). (C.l) 
Under this transformation, N, Aj and gij are found to transform as 

= A, - diS^Xj - eXi - die (N^ + A^) - ^,,■e^ 
9ij = 9ij - die^9kj - dj^^Qik - die Xj - dje Aj. (C.2) 
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Furthermore, F*-;,, the affine connection defined by Qij, transforms under the diffeomor- 
phism (|C.1|) as 

rj. = r),-d, d,e^ + F- dme' - Tl, d,e^ - F}^a,e™ + 5F}„ (C.3) 

with 

~5V% = -VV. Vfce - S.eVfc Y - a,.eV, Y - Ng'\d,e Kiu + due Kij - die K^u). (C.4) 

Note that SV^j^ does not contain e*. From these relations, it is straightforward to verify 
that the extrinsic curvature transforms as 

K[^=K,^-d.,e'Ki^-due'Kji 

+ NViVje + d^e {djN - X'Kji) + {d^N - X'Kij). (C.5) 

We can also show that the Riemann curvature R^jki transforms under ( |C.1| ) as 

r)'i _ p* j_ f) J P"^ Pi T^i f) ^"^ P* f) ^"^ P* 

^ jkl — ^ jkl "T ^me JTL jj^i Ujt JTL ^j^i Cfee JTL jj^i Ult JTL jf^^ 

-due % + die tl^ + Vfe5Fi^. - Vi~5Ti^. (C.6) 



As argued in §3, we focus on the diffeomorphism that obeys the condition ( p.6| ). This 
is equivalent to the following relation in an infinitesimal form: 

9ie(r = ro) = 0. (C.7) 



Therefore, we find that the boundary action ( |3.7|) is invariant under this diffeomorphism. 
We remark that in the above, we have discarded boundary terms of the form 

S',= I d''x^{K'^U^ + Kg'^U^), (C.8) 

although these are allowed by the diffeomorphism.^ The reason is that if there were such 
boundary terms, they would require us to further introduce an extra boundary condition, 
since 

5S[ = [ d^x^ [■■ + 5k,,Pi\gku Kki)] . (C.9) 



^^By definition, the {d + l)-diniensional scalar curvature i? is a scalar. It thus follows from (B^) that 
Lij{r = rQ) transforms as a tensor under the diffeomorphism with 
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D Coefficients in Eq. ( ^:t8l) 



We have the following values for the coefficients in Eq. (|3.18|) : 



Ai 
A, 
A, 



26 -X2 



A,, 



46 + 8c - 3x5 



A, 



b + Xi 



2(6 + 4c)' 2(6 + 4c) ' 6 + 4c' 

4a6 — 16ac + hxi + 4cxi — Aax2 + 26^ — hx2 
2(6 + 4c)(4rfa+ (d + l)6 + 4c) ' 
4a6 — 16ac + 26^ + 6x4 + 4cx4 — 12ax5 — 36x5 
2(6 + 4c)(4rfa+ (rf+ l)6 + 4c) ' 
4a6 — 16ac — 36x3 — 12cx^ + 80x4 + 26^ + 26x4 
2(6 + 4c)(4da+ (ci + l)6 + 4c) ' 



1 



4(6 + 4c)(4rfa+ (d+ l)6 + 4c) 



X 



46^ + 4((i + l)a6^ + 4ax2^ - 46^X2 + 6x2^ + 64ac^ - 806x2 



+16((i — 2)a6c — Adcxi' — dbxi^ + 46^Xi + 166cxi — 8cxiX2 



+32acx2 — 26x1X2 + Sdabxi + 32dacxi 



166c + 46x2 — 
4(6 + 4c) 



5, 



c, 



(D.l) 



(D.2) 



1 



4(6 + 4c)(4c/a+ (c/+ l)6 + 4c) 

X 86^ — 8a6x2 — 16{d + l)6cxi — 64c^Xi — 32dacxi 

—Adb^Xi + Sdabxi — 46^X2 + 32acx2 + 8 dabx^ — 2Aabx^ 

+24ax2X5 + 66x2X5 — 126^X5 + 32((i — 2)a6c — 2dbxiXi 

+9>{d + l)a6^ + I66CX4 + 46^X4 — 26x2X4 — 66x1X5 

— 8CX2X4 — SdcxiXi + 32(iacx4 — 24cxiX5 + 96acx5 + 128ac^ 
1 

4(6 + 4c)(4c/a+ (c/+ l)6 + 4c) 

X - 166^c + 86cx2 + 64c^xi + 32dacxi - 32(rf + 2)a6c 
— 8(7(i + 5)a6^ — {d — 3)6x2X4 — A{d — 4)0x2X4 



(D.3) 
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+S{d - 2)06x4 + 2(d - 3)6^x4 + 3{d - 1)6x2X3 - 6{d - 1)6^X3 
-A{d + 3)b^ + 32acx2 + 2A{d + 1)06x3 + 16(rf + l)6cxi 
+64(ia^X2 - I2CX2X3 + 2{d + 3)6^X2 + Sdabxi + Adb'^xi 
—ddbxiXs + 246cx3 — 46x1X4 + 12^0x2X3 + 860x4 + 96dacxs 
— 4CX2X4 — I6CX1X4 + 64acx4 — 24dcxiX3 — 128da^6 — 128ac^j , 

326c + 66x5 — 3x2X5 + 64c^ — 8cx2 
C3 = 



a 



2(6 + 4c) 

—86c + 26x4 — 26x1 — X2X4 — 8cxi + 4cx2 



D3 



6 + 4c 

C5 = -12c-2x2, (D.4) 
86c - 9x5^ - 48cx5 - 32c^ 



4(6 + 4c) 

—46x5 ~ 166c — I6CX4 — 6x4X5 + 8CX5 



2(6 + 4c) 
1 



4(6 + 4c)(4c^a+ (c^+ l)6 + 4c) 



X 



66x4X5 — 64:C^X4 + 96acx5 — 16{d + l)6cx4 — 32dacx4 + 128c^ 

-4ci6^X4 - 24CX4X5 + 32(0? + 2)6c^ - dbxi^ - Adcx^ + 8(ci + l)6^c 

+46^ + 64(2d + l)ac^ + 4(d + l)a6^ + 16(3d - 2)a6c - 24a6x5 

— 126^X5 — Sdabx4 + 96x5^ + 36ax5^ 
1 



4(6 + 4c)(4da+ (d+ l)6 + 4c) 



X 



86'^ — 32cx4^ + 480x4X5 — 16dax4^ + 24a6x5 + 126^X5 + 126x4X5 
—2Adabxz — 96ciacx3 + 64c^X4 — 96acx5 — 192c^X3 
-72CX3X5 + 32{d + 2)a6c - 48(d + l)6cx3 - 8(3d + 2)06x4 
+16(c^ - l)6cx4 + 32(ci + 2)00x4 + 16(rf + 1)6^0 - 4(rf + 2)6x4^ 
-A{d + 4)6^X4 - 8((i + l)a6^ - 128ac^ - 6^6x3X4 
-186x3X5 - 24dcx3X4 - 12d6^X3 + 646c^] , 
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4(6 + 4c)(4rfa+ l)6 + 4c) 

X IQax/^ + 64c^X3 — 8dabx3 — 32dacx3 — 126x3X4 — 480x3X4, 
+A{d - 2)6^x3 - 64acx4 + 86^X4 + 46x4^ + 46^ + QAac^ - 9dbxi^ 
-36c/cx3^ + A{d + l)a6^ + 16(rf - 2)ahc + 1606x4 + 16(rf - 1)60x3 



(D.5) 



El = 46 + 2xi - X2, 
E2 = -2b + 8c -2x2, 
E^ = -8c - X2, 
E4 = 26-2x1. 



(D.6) 
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